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Sommario. In questo lavoro abbiamo studiato un modello di omogeneizzazione per una 
classe di operatori tipo Kolmogorov-Fokker-Planck in forma di divergenza. Le proprietà 
di questa classe di ‘operatori sono state oggetto di una ricerca-sistematica solo negli ultimi 
anni, a partire dal lavoro di Lanconelli e Polidoro [15]. In questo lavoro, in particolare, 
vengono studiate le proprietà fondamentali della classe di operatori di evoluzione considerata, 
nel caso in cui i coefficienti soddisfino alcune ipotesi: la classe di operatori considerata è 
di tipo Hérmander e risulta invariante rispetto ad un particolare gruppo di traslazioni e 
rispetto ad una certa famiglia di dilatazioni. L’ esistenza di tale gruppo omogeneo è |’ 
ingrediente basilare per sviluppare uno schema di omogeneizzazione. Infatti la teoria esistente 
sull’ omogeneizzazione si basa essenzialmente sull’ azione delle usuali traslazioni e dilatazioni 
di tipo euclideo in ®". In questo lavoro, invece, abbiamo sviluppato una teoria analoga 
basata, però, sulla esistenza di altri gruppi ed adattata quindi ad operatori invarianti rispetto 
a tali gruppi come in Biroli-Mosco-Tchou (si vedano [28] e [29]) che hanno considerato un 
problema analogo, nel caso stazionario, per operatori ellittici degeneri associati al gruppo 
omogeneo di Heisenberg. In sostanza la costruzione di un particolare pavage dello spazio 
utilizzando il gruppo di traslazioni e di omotetie associate con la classe di operatori tipo 
Kolmogorov-Fokker-Planck considerata, permetterà, in analogia con la teoria classica dell’ 
omogeneizzazione, la formulazione di un problema di omogeneizzazione periodica. Inoltre 
attraverso un procedimento di regolarizzazione, definito come regolarizzazione subellittica, 
riusciamo a lavorare nel contesto delle soluzioni deboli analogamente al caso classico. Al 
contrario di questo ultimo, non applichiamo il metodo standard dell’ energia per provare i 
voluti risultati di convergenza nel procedimento di omogeneizzazione, ma ci basiamo su un 
metodo alternativo, noto con il nome di sviluppo aggiunto, il cui ingrediente fondamentale è la 
costruzione di una classe di funzioni test grazie ad un particolare metodo di regolarizzazione 
dei coefficienti dell’ operatore studiato. Infine otteniamo il voluto risultato di convergenza ma 
a meno di una sottosuccessione, a causa dell’ assenza di risultati di unicità per la soluzione 
debole del problema omogeneizzato. 
Abstract. In this paper we study an homogenization model for a class of Kolmogorov- 
Fokker-Planck operators in divergence form. The features of this class of operators were 
subject of a systematic research only in the last years, starting from Lanconelli and Polidoro 
work [15]. In that paper, in particular, the essential properties of the considered class of 
evolution operators are studied, when the coefficients satisfy some conditions: indeed the 
class of operators is of Hérmander type and it is invariant with rispect to suitable groups of 
translations and dilations. The existence of such an homogenous group is the main ingredient 
in order to develop an homogenization scheme. In fact, as we will show in the first part of 
this paper, the existing homogenization theory refers to the action of the usual translation 
and dilation of euclidean type on R!. On the other hand, in this work, we have developped a 
similar theory which requires the existence. of other groups and suitable for operators invariant 
for such a groups, following the Biroli-Mosco-Tchou ideas (see [28] and {29]) who considered 
a similar problem, in the stationary case, for degenerate elliptic operators associated with 
the homogeneous Heisenberg group. Substancially the costruction of a particular pavage of 
the space using translations and dilations group associated with the class of Kolmogorov- 
Fokker-Planck operators, will allow, as in the classical homogenization theory, the definition 
of periodic homogenization problem. Furthermore using a procedure of regularization, defined 
as subelliptic regularization, we are enable to consider weak solutions, as in the classical case. 
Despite the latter case, we don’ t use the standard method of energy to proof the desired 
convergence results, but we will use an alternative method, known as adjoint erpansion, 
which the most important feature is the contruction of test functions class using a peculiar 
regularization method of the coefficients of studied operator. Finally we obtain the expected 
‘ convergence result but up to a subsequence, because the absence of unicity results for weak 
solution of homogenized problem. 
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1 Introduzione 


I materiali compositi (fibrati, stratificati, porosi,...) giocano un ruolo molto importante in 
vari campi della Meccanica, della Fisica, della Chimica e della Ingegneria. In termini molto 
generali, un mezzo composito è un mezzo altamente non omogeneo e non isotropo, con delle 
eterogeneità al suo interno che ne alterano da luogo a luogo le proprietà fisiche rilevanti. 
Quindi i parametri fisici (come la conduttività, i coefficienti di elasticità...) di questi mate- 
riali sono tipicamente discontinui ed oscillano tra i differenti valori che caratterizzano ognuno 
dei componenti. Quando questi componenti sono finemente miscelati, questi parametri oscil- 
lano molto rapidamente e la struttura microscopica diventa abbastanza complicata. Pertanto 
questi materiali sono fortemente eterogenei, se osservati ad un livello microscopico, ma esi- 
biscono il comportamento tipico di un materiale omogeneo ideale, quando sono osservati su 
scala macroscopica. a 

Allora, possiamo pensare di ottenere una buona approssimazione del comportamento 
macroscopico di questo materiale eterogeneo, mandando a zero il parametro €,, che de- 
scrive la finezza della struttura microscopica nelle equazioni che intervengono in fenomeni 
come, ad esempio, la conduzione del calore o l’elasticità. Questo è lo scopo della teoria dell’ - 
omogeneizzazione, ossia descrivere tali processi limite, quando ex tende a zero. . 

Più precisamente, l’omogeneizzazione si occupa dell’ analisi asintotica di equazioni alle 
derivate parziali governanti alcuni processi fisici che intervengono in materiali eterogenei con 
una struttura periodica, quando la lunghezza caratteristica del periodo €, tende a zero. 

Un buon modello per lo studio del comportamento fisico di un materiale eterogeneo con 
una struttura periodica molto fine, come in elettrostatica, in magnetostatica, o quando si 
considera la parte stazionaria della diffusione del calore, è dato dal seguente 


-div(o(£)Du) = finQ, ti 
unloo = 0su09; 


qui £ è un insieme limitato di PR" che consideriamo come una parte del materiale eterogeneo 
ed € è il periodo della struttura in tutte le direzioni. 

La funzione u, può essere interpretata, in base al campo di applicazione, come il poten- 
ziale elettrico, oppure come il potenziale magnetico, o come la temperatura. La matrice 
dei coefficienti a(£) = (a;;(£)) è costituita da funzioni €y-periodiche e descrive le proprietà 
fisiche dei differenti materiali costituenti il corpo preso in esame (essi sono i coefficienti dielet- 
trici, ad esempio, o la permeabilità magnetica oppure i coefficienti di conduttività termica). 
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La funzione f è un termine di sorgente assegnato. Un semplice esempio, per chiarire le idee, 
è fornito dal miscuglio fine di due materiali della stessa natura ma di caratteristiche quanti- 
tativamente diverse, come due materiali conduttori di conducibilità diverse c, e c.. Un mezzo 
di questo tipo può essere descritto da un tensore di coefficienti caratteristici a;;(r) = 4;;(7) 
i,.j = 1,...,n, discontinuo nella variabile x € R” , con a;;(r) = ci nella parte del corpo 
occupata dal materiale 1 e con a;;(r) = c. nella regione occupata dal materiale 2. Quindi, 
ritornando al nostro problema, se il periodo della struttura è molto piccolo, una approssi- 
mazione numerica diretta della soluzione del problema (1) può essere molto difficile se non 
impossibile. 

Pertanto l’omogeneizzazione ci fornisce un modo alternativo di approssimare tali soluzioni ux 
tramite l'utilizzo di una funzione u che risolve il problema corrispondente ad un materiale 
omogeneizzato 


—div(bDu) fin9Q, 
[ ulsa = 05su 09, (2) 


dove b è una matrice costante (ricordiamo che per un materiale omogeneo le proprietà fisiche 
non dipendono da x). La matrice omogeneizzata b rappresenta i parametri fisici di un corpo 
omogeneo, il cui comportamento è equivalente, da un punto di vista macroscopico, a quello 
del materiale con microstruttura periodica assegnata (descritto tramite la (1)). Quindi la idea 
è quella di sostituire, in una prima approssimazione, il comportamento del materiale reale 
non omogeneo con struttura e, - periodica con quello ideale omogeneo e dunque sostituire le 
soluzioni u, con la soluzione u relativa al problema (2) più idonea ad essere trattata da un 
punto di vista numerico perché soluzione di una equazione a coefficienti costanti. 

Le principali referenze per la teoria dell’ omogeneizzazione di strutture periodiche sono 
i libri di Bensoussan-Lions-Papanicolaou [1], Sanchez Palencia [2], Lions {3], Bakhvalov- 
Panasenko [4], e Jikov-Kozlov-Oleinik [5]. Per quanto riguarda i principali risultati astratti 
di analisi funzionale glie verranno sistematicamente utilizzati nel seguito (convergenza de- 
bole, spazi di Sobolev e loro proprietà, lemma di Lax-Milgram, disuguaglianza di Poincaré) 
ci riferiamo al testo di H. Brezis [6]. 


2 Definizione del problema e costruzione dell’ operatore omo- 


geneizzato 
Siano z1,...,2n n vettori linearmente indipendenti di R", e sia P il parallelotopo con lati 
Zix-++yZn, CIOÈ 


P= ftiz:+---+tnzn :0<t<1iî=1,...,} 


Diciamo che una funzione  : R" — R" è P-periodica se (7) = «(x + 2z;) per ogni x E R" e 
per ogni i = 1,...,n. In tal caso noi diciamo che P è una cella di periodicità della funzione 
«. Per semplicità (ma senza perdita di generalità) noi assumeremo sin da ora che la cella di 
periodicità P sia il cubo unitario Y =]0, 1[°. Quindi se {ei,...,en} denota la base canonica 
di R" allora 4 : R" + R" è Y-periodica se p(r) = g(r + e;) per ogni © € IR" e per ogni 
= ssh n 

Denotiamo con My l’insieme di tutte le funzioni a : R"° + M?°*", con a(z) = (a;;(7)) per 
x € R°, soddisfacenti le seguenti proprietà: 


(1.2.1) a;; è Y-periodica in R" per ogni î,j=1,...,; 
(1.2.2) a;; € L°(R") per ogni ig=l, uf 


(1.2.3) esiste una costante a > 0 tale che (a(x)€,€) = Vi;=i ij (2)6;G > a|€?| per quasi ogni 
i t ER" e per ogni £ E R”. 
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Definiamo a’ : R" — M"*" come x 
al (e)= 2), (3) 


dove €, è una successione di numeri reali positivi convergenti a zero. Si osservi, inoltre, che 
le funzioni al. sono €xY-periodiche in IR". Sia 2 un dominio limitato ed aperto di R" (lo 
considereremo come un pezzo del materiale eterogeneo in esame). Per un fissato EL > 0, 
consideriamo il seguente problema di Dirichlet ai valori al contorno su Q 
-div(a*(2)Du) = finQ, 4 
urlan = 0su09, (4) 
dove f è una funzione assegnata su 9. Assumiamo che f € H°!?(9). La formulazione 
variazionale (debole) del problema (4) diventa: trovare u, € H3'(9) tale che 
Sala" (2)Dux, Dv)dz = (f,v) per ogni v € H;° (0) 
[ un € Hi" (0) Là 
(si osservi che questa presentazione non richiede ipotesi di regolarità sui coefficienti dij; co- 
munque, come vedremo, noi otterremo stime a priori su u, che sono indipendenti da £, e che 
non si basano sulla regolarità dei coefficienti). Grazie al Lemma di Lax-Milgram si ottiene |’ 
esistenza e l’ unicità della soluzione del problema (5). Dunque da un, punto di vista teorico 
il problema è completamente risolto. ; : 
Per studiare il comportamento asintotico delle soluzioni un del problema (4), dove supponi- 
amo f € L?(Q), una tecnica efficiente consiste nell’ applicare lo sviluppo asintotico utilizzando 
scale multiple (ossia, variabili lente e veloci). Rimandiamo per questa tecnica ai testi citati, 
in particolar modo a quello di Bensoussan-Lions-Papanicolaou [1] ed a quello di Sanchez- 


Palencia [2]. Più precisamente si suppone, in modo del tutto euristico, che u ammetta uno 
sviluppo in due scale della forma 


un(2) = wo(2, 2) +eyu (2,2) +efu(1,Z)+..., (6) 
Eh Eh Eh 


dove le funzioni u;(2,y) sono Y-periodiche in y per ogni x € Q. Questo significa che noi 
postuliamo l’ esistenza di funzioni regolari u;(r,y) definite in Q x R", Y-periodiche in y ed 
indipendenti da e, in modo tale che il lato destro della (6), per y= (=); sia uno sviluppo 
asintotico di u,(z). 

Notiamo che le due variabili x e (#) prendono in considerazione le due scale del fenomeno 
della omogeneizzazione; la variabile £ è la variabile macroscopica, mentre la variabile (£) 
si riferisce alla geometria microscopica. Il metodo che stiamo andando a sviluppare si riv- 
elerà molto utile per ottenere le giuste risposte sul comportamento al limite delle soluzioni 
del problema (4). Facciamo osservare, per ultimo, che è tecnicamente complicato tenere in 
considerazione le condizioni al contorno quando cerchiamo la funzione u, nellà forma (6), e 
questo rappresenta al momento una sorgente di serie difficoltà tecniche per la giustificazione 
del metodo. Per maggiori dettagli su queste difficoltà si rimanda all’ articolo di Sanchez- 
Palencia [7]. Il metodo, noridimeno, fornisce le risposte giuste poiché in questo particolare 
tipo di problemi molto spesso le condizioni al contorno sono in qualche modo irrilevanti. L’ 
idea del metodo è (semplicemente) quella di inserire la (6) nella (4) ed identificare le potenze 
successive di ey. Per presentare questi calcoli in una forma molto semplice, data una funzione 
regolare $(7,y) di due variabili, definiamo la funzione ®, (x) tramite la 


D:(2)= (2,2) 
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ed osserviamo che EEA BI) 10% 
h/-y_ 09 100, © 
da p (Fa, le Eh dy) P E 


Allora, si può scrivere 


AS, = [(ep2A0 + 67141 +84) 9](2, —), (7) 
dove 
o 0 le) 
Ao= L By CH 0) da) (8) 
2 0 d 2. 0 d 
A =- pa dal 0) = p> By (0907) (9) 
A:=- D 7204075) (10) 


Utilizzando la (6) e la (7), I equazione (4) è soddisfatta, sotto l’ ipotesi che a e le funzioni 
u;(x,y) siano sufficientemente regolari, se 


(Asuo)(r,=)=0in9, (11) 
Eh 
(Aot + Aruo)(r,) =0inQ, (12) 
h 
(Aouz + Art + Azto)(£, >) = finQ; (13) 


cerchiamo adesso di capire come dalle (11), (12) e (13), si possa costruire |’ operatore omo- 
geneizzato. Determiniamo le funzioni w; risolvendo i seguenti problemi nei quali consideriamo 
T come un parametro fissato: 


(Aouo)(2,:)=0inY, 
{ ir Y-periodica ; (14) 
(Ao )(7,:) = —(A1vo)(2,-)in Y, 
{ u(z, .) Y-periodica ; i (15) 


{ (Aouz)(2,:) = f() — (Aru + Azwo)(2,-) in Y, (16) 


uz(,-) Y-periodica . 


Consideriamo questi problemi nel contesto delle soluzioni deboli. Cominciamo col ricordare 
un risultato di esistenza di soluzioni per un problema ai valori al contorno sul cubo unitario. 
Denotiamo con H};?(Y) il sottoinsieme delle funzioni u di H*:?(Y) che hanno la stessa traccia 
su facce opposte di Y. Inoltre indichiamo con HEY) il sottoinsieme di H!:?(Y) di tutte le 
funzioni con valore medio nullo e che hanno stessa traccia su facce opposte di Y. 


Lemma 1 Sia F € (H}:?(Y))". Allora il problema 


per 


Sy(a(y)Dy9, Dyb)dy = (F, 6) per ogni w € Hxi:(Y) , 1”) 
de Hg°(Y) 


ammette una unica soluzione se e solo se 


fF@dy=o (18) 
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Per la prova di questo lemma rimandiamo a [1], pag. 13. Sfruttiamo questo lemma per 
determinare le soluzioni dei problemi (14), (15) e (16). 
Step 1: Studio del problema (14). 

Utilizzando la formula di Green e ricordando le ipotesi di periodicità, si prova facilmente 
che il problema (14) può essere formulato in modo debole come segue: trovare Uo(7,-) E 
H}:(Y) tale che 


per 
VICOLI Dyw)dy= 0 per ogni w € Haie(Y). (19) 


Ricordiamo che x gioca il ruolo di un parametro e che Dyuo è il gradiente nella variabile 
periodica y. Allora, grazie al Lemma (1) possiamo concludere che uo(£,-) è determinata da 
(19) a meno di una costante. Prendendo v = wo(£, +) nella (19) ed usando la condizione di 
ellitticità soddisfatta da a segue immediatamente che Uo(£, -) =costante, cioè 


Uo(2, +) = vo(2). (20) 


Step 2: Studio del problema (15). 
Utilizzando la (20), il problema (15) si riduce al seguente 


(Aov)(2, )= Dig (ap ;(1)) fe (2) int, (21) 

u1(, -) Y-periodica. 
Questo è ancora un problema in y, dove x è un parametro. È possibile rappresentare w; in 
una forma molto semplice grazie alla separazione delle variabili nel lato destro della (21). 
Infatti, grazie alla formula di Green ed alle ipotesi di periodicità, la formulazione debole del 
problema (21) diventa: trovare w1(7, +) € H1:? (Y) tale che 


per 
= duo - x71,2 
feta, D,pray=- DD / asus ENO. (22 
i,j=1 
Allora consideriamo per j = 1,...,n il problema 
Sy (a(Y) Dy wi (y), DyW(y)) = — fr (a(y)e;, Dyb(y))dy vwe Hy°(Y), 9 
w e Hy°(Y). (23) 


Si osservi che l’ applicazione Fi : Hi°(Y)+R definita da 
vò FI(4)=— f (a(%)e;, D,)dy 


è un funzionale lineare e continuo su BE). Pertanto essendo Jo (Gr) dy = 0, in virtù del 
Lemma (1), esiste ed è unica la soluzione del problema (23) per ogni j = 1,...,n. Allora la 
soluzione generale del problema (22) diventa 


(2,9) = Dei a, (9), (24) 


dove %, è una costante additiva (dipendente, quindi, dal parametro 2). 
Step 3: Studio del problema (16). 

Passiamo ora al problema (16) dove consideriamo u, come l’ incognita, essendo x un 
parametro. Ancora una volta, in virtù del Lemma (1), u2 esiste se e solo se 


(Ax, + Asto)dy = i f(2)dy= |Y]f, (25) 
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dove |Y| indica la misura di Y. La condizione (25) è proprio equazione omogeneizzata che 
stiamo cercando. Si noti che 


di lo) du, 
Apiadpo= SÈ J (Mera 
J e Zia 900 pu, i 
ed utilizzando la (24) si ha che 
= 0 2 Qwi duo 
Agg E J at mago 
f 1 dy z Fra 2000) du, VU Fa. 
per cui la (25) si esprime come 
ui 0? uo 
- bi ——-=finQ, (26) 
z 19x;07; 
dove dwi(y) 
1 o wi (y 
bi; = 7J i ly —_)d 27 
i [Y] |, (Gi; (0) + 2a 0 (0) dY ) Y ( ) 


L’ equazione (26) è l’ equazione omogeneizzata mentre i coefficienti b;; sono i coefficienti omo- 
geneizzati. Infine per ottenere un problema ben posto per uo bisognerà soltanto imporre una 
condizione al contorno per uo. Dalla (4) e dalla (6) si ricava che ulsn = 0 su 09. Si osservi 
che questa relazione è formale ma, come vedremo, è rigorosamente giustificabile. 


Osservazione. Possiamo riassumere le precedenti considerazioni e concludere che se si ipo- 
tizza uno sviluppo della forma (6), il primo termine è determinato come la soluzione dell’ 
equazione (26) con la condizione al contorno volan = 0 su 02. Inoltre il procedimento formale 
(che verrà giustificato nel prossimo paragrafo) per calcolare i coefficienti omogeneizzati è il 
seguente: 


i) si risolve la (23) sulla cella unitaria Y,perg=1aaa5 
ii) i coefficienti omogeneizzati bi; sono dati dalla (27). 


Nel prossimo paragrafo mostreremo che la successione (un) converge, in una appropriata 
topologia, alla funzione vo sopra assegnata. 


3 Omogeneizzazione 


Data a € My, consideriamo il seguente problema di Dirichlet ai valori al contorno sul sot- 
toinsieme aperto e limitato Q di R": ° 
{ -div(a(£)Du) = finQ, 


u € Hl*(9), (28) 


dove f € H-!:?(Q) ed €, è una successione di numeri reali positivi convergenti a 0. Conside- 
riamo inoltre il seguente problema omogeneizzato 


-div(bDu) = finQ, 
{ uo € HI (9). (20) 
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in cui la matrice b = (b;;) è definita tramite la 


250) gy, (30) 


(i) = fl (i;(0) +Y0 a(0) 
x k=l 
dove wi (y) è I’ unica soluzione del problema locale 


Je (a(y)(e; + Dywi (y)), Dyv(y))dy = 0 per ogni v € Hy°(Y) E 31 
wi e Hy°(Y). (31) 


Allora è valido il seguente 


Teorema 2 Sia a € Ma ed (ex) una successione di numeri reali positivi convergenti a 0. 
Siano un ed uo rispettivamente le soluzioni dei problemi (28) e (29). Allora, 


Un + uo in H3'"(Q)-debole , (32) 
a(—)Du + bDuo în L°(Q;R")-debole. |. (33) 
h a 


Per la prova di questo teorema si rimanda al testo di Bensoussan-Lions-Papanicolaou [1] 
(Cap.1, Teorema 3.1). Vorrei invece sottolineare alcune caratteristiche di tale dimostrazione. 
Questa si basa su un metodo molto potente e generale conosciuto col riome di metodo dell’ 
energia introdotto da Tartar. Il suo nome, dovuto al fatto che rappresenta la prova diversa di 
un importante teorema di convergenza di De Giorgi e Spagnolo (vedi [8],[9]), non è realmente 
adeguato, poiché il suo principale ingrediente è una accurata scelta di funzioni test. Più 
precisamente lo scopo fondamentale di questo metodo è di passare al limite nella formulazione 
variazionale del problema (4): 


/, a'(2)Vun(2) - Vé(z)dz = J, f(2)$(x)dz per ogni $ € H1?(9) . (34) 


Per una data funzione test $ non si può passare al limite nella (34), quando €, tende a 
zero, poiché il lato sinistro coinvolge il prodotto di due successioni debolmente convergenti (si 
può infatti mostrare che a’ converge, in L?(G)-debole, alla funzione M(a) = a] Sa(y)dy). 
Allora l’ idea principale del metodo è quella di sostituire la funzione test fissa $ con una 
successione di funzioni $, accuratamente costruita che permette di passare al limite grazie 
al fenomeno della compattezza per compensazione (vedi [10] e [11] per questa nozione). La 
giusta successione di funzioni test è 


E) = 0 +a DT E), (65) 


dove $ E C5°(0) e wi è la soluzione del problema aggiunto locale (cioè il problema (31) in 
cui la matrice a è sostituita dalla sua trasposta che indichiamo con a*). Integrando per parti 
nella (34) ed utilizzando la equazione (23) sulla cella unitaria riusciamo a passare al limite 
nella (34) ottenendo la formulazione variazionale del problema omogeneizzato (29). In questo 
modo la convergenza del processo di omogeneizzazione è rigorosamente provata. 

Sottolineo, ancora una volta, |’ importanza rappresentata da una scelta accurata di fun- 
zioni test, perché, come vedremo successivamente, questa scelta sarà di fondamentale impor- 
tanza per provare i risultati di convergenza nel procedimento di omogeneizzazione della classe 
di operatori tipo Kolmogorov-Fokker-Planck studiati. Infine è importante ricordare anche il 
seguente risultato, provato sempre in [1], grazie al quale il problema (29) risulta ben posto: 


96 


Proposizione 3 Sia a : R" + M"*" una funzione di My. Sia b la matrice costante i cui 
coefficienti sono definiti dalla (1.34). Allora b soddisfa la stessa condizione di ellitticità di a, 
cioè, 


(68,6) = DO bi;6;6; > alé|? per ogni € ER". (36) 


i.j=1 


Adesso possiamo passare ad analizzare il nostro problema, ovvero lo studio di un modello di 
omogeneizzazione e di regolarizzazione per una classe di operatori tipo Kolmogorov-Fokker- 
Planck. Per la prova dei principali risultati si rimanda al lavoro [12]. 


4 Preliminari 


Cominciamo col richiamare alcune proprietà della classe di operatori tipo Kolmogorov-Fokker- 
Planck che noi considereremo, dettagliatamente studiate nel lavoro già citato di Lanconelli e 
Polidoro [15]. Attualmente, lo studio di questa classe di operatori ha aperto un vasto campo 
di ricerca ricco di profonde e delicate problematiche. Un elenco non esaustivo di contributi 
importanti può essere dato dai lavori [13], [14], [16], [17], [18], [19], [20], [21], (22), [23], [24], 
[25], [26] e dai lavori citati nelle relative bibliografie. In particolare, nel lavoro di Lanconelli e 
Polidoro [15], non solo viene individuato un gruppo di traslazioni per cui la classe di operatori 
risulta invariante ma vengono sottolineate anche le condizioni da imporre (vedi ipotesi (H2) 
sotto) affinché esista un gruppo di dilatazioni rispetto al quale l’ operatore sia omogeneo di 
grado 2. 
Focalizziamo quindi la nostra attenzione sulla classe di operatori di evoluzione del seguente 
tipo: A 
L= > 9x;(0;;(2)0:,) + (©, BD)-d (37) 
- dgj=1 
dove 2 = (x,t) € RVH!, D= (9:,,---19:y) ed (-,:) è I usuale prodotto interno di RN. Sui 
coefficienti di L faremo le seguenti ipotesi: 
(Hi) (a;;(2)) = (a;:(2)) € L°(RN+!) per ogni î,j= 1,...,9 ed esiste p > 0 tale che 
q 


De < Das <n DE (38) 


i.j=1 i= 


rlW 


per ogni z € RN+! e per ogni (£1,...,£3) E R®. 


(H.) B è una matrice costante N x N avente la seguente forma: 


0 B;, 0 0 
0 0 BB. 0 
IE $i (39) 
0 0 0 ... B 
o 0 0 
dove per ogni i = 1,...,r B; è una matrice p;-1 X p; con rango pi € A = Po > pi>...> 


Pr Po+pi+::+p= N. 


Come mostrato în [15], Proposizione A.1, le condizioni (H1) ed (H2) assicurano che l° oper- 
atore a coefficienti costanti, cioè dove a;; = a; sono termini costanti per ogni tg = 1:+434 
ancora verificanti la (38) 


Lo = SO 8,,(0159,,) + (3, BD) - 2 (40) 


i.j=1 
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è ipoellittico. Infatti, posto 
Y=(x,BD)-d, 


uno dei risultati principali di [15] consiste nel provare che la condizione (HA) è equivalente 
alla seguente: 


(H3) rango Lie(0,,,...,0x,,Y)(2) = N+1Vz € RNH dove Lie(9:,,-..,0z,,Y) indica l” 
algebra di Lie generata dai campi d,,...,0:,,Y. 


Come è ben noto (vedi Hòrmander [27]), la condizione (H;) implica la ipoelliticità dell’ oper- 
atore (40). Facciamo presente che d’ ora in avanti ci rifereremo all’ operatore con coefficienti 
costanti tramite la notazione Lo in consistenza con la (40). 

Sempre in [15] si verifica che 


(2,0) 0 (4,7) = (y+ E(7)2,t+7) per ogni (2,8), (y,7) € RN#! (41) 


dove E(r) = erp(-r.BT), allora (RN+, o) è un gruppo il cui elemento neutro è (0, 0); 1’ inverso 
di un elemento (7,t) € RNH! è (2,t)-1= (-E(-t)x, -t). Ogni operatore Lo è invariante per 
le traslazioni a sinistra di R#, cioè definite le traslazioni a sinistra come segue 


li :RNY a RNH, 1,.(2)=2z'oz (42) 
dove 2' € RN+! è fissato, si verifica immediatamente che 
Lo(u(2'0 z)) = (Lo(u))(z"0 z) (43) 


per ogni u € C8°(RN#1) e per ogni 2 € RNH con ' fissato. 

Inoltre, sempre in [15], si mostra che la particolare forma della matrice B indicata nella 
(H2) assicura l’ esistenza di un gruppo di dilatazioni su RN+, indicato con (D(A))\.0 rispetto 
al quale Lo è omogeneo di grado 2 cioè: 


Lo o D(A) = A°D(A) o Lo. (44) 
La dilatazione D(A), grazie all’ ipotesi (H2), è la seguente 
D(A) :RN+ 4 pH, D(A)(2,t)=(Az;,.. « A°NEN, A°N41t) (45) 


dove a, = ---= ap, =], &po+41 °°° = Qpo+pi = 3, --.3 @po+--4p,_141 == N =2F+1 ed 
@n+1 = 2. Possiamo quindi scrivere 


D(A) = diag(Mpo, A3Ip,xAIpjr..., ASSELE AR) (46) 


dove I, è la matrice identità k x k. In forma più esplicita possiamo riesprimere la (44) nel 
modo seguente: 
Lo(u(D(A)z) = A?(Lou)(D(A)z). (47) 


Si tenga presente che esistono due proprietà molto importanti che tollegano le dilatazioni 
D(A) e la legge di gruppo (R+!,0) associata ad Lo: 


D(A)(2'0 2) = (D(A)2') 0 (D(A)2) per ogni 2,2" € RN, (48) 


D(A)(27!) = (D(A)2)7! per ogni 2 e RN+!. (49) 


Per finire, ricordiamo la seguente definizione di una norma in RN+!, omogenea di grado 1 
rispetto al gruppo di dilatazioni (D(A))x30, introdotta nel lavoro di M. Manfredini [21]. 
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Definizione 4 Siano 01,...,@n le costanti per cui la matrice D(A) si scrive nella forma 
(46). Se 2 € RN+!\{0} poniamo || 2 ||= p se p è l’ unica soluzione positiva dell’ equazione 


e||2]|=0sez=0. 

Inoltre, ancora in [21], è verificata la seguente 

Proposizione 5 L' applicazione z+|| 2 || ha le seguenti proprietà: 

Ni) || D(A)z ||= A || 2 || per ogni 2 € RN! e per ogni A > 0; 

N») L’ insieme {z/|| z ||= 1} è la sfera euclidea Yy41 = {(2,t)/lc}+|i]}=1}. 


N3) Esiste una costante co = co(B) > 0 (dipende solo dalla matrice B) tale che per ogni 
z,( E RN+ si ha 
zo 6 IS coll 2 {+1 61: 


N.) Esiste una costante c = c(B) > 0 tale che 


1 
22 <N27* {<a Il 21 per ogniz e R*4; 
3 È 


Ns) Per ogni compatto K di RN+! esiste c, = c2(B, K)> 0 tale che 
|l:- | <c2|| oz] per ogni 2,6 € K con || Coz|<1, 
dove | -| indica la norma euclidea in RN. 
Sempre seguendo [21], introduciamo una quasidistanza in RN+! relativa alla norma || ||: 
Definizione 6 Per ogni 2,( € RN#! poniamo 
d(2,0)=|| oz. 

d è una quasidistanza (Cfr. [21]), nel senso che 

d(2,6) 2 0, d(z2,60)=0seesolosez=(C, 
ed esiste c > 0 tale che 

d(2,6) < c(d(2,29) +d(2°,0)) Vz,24,6 € RP. 
Denoteremo con B(z,r) l’ insieme 
B(z,r) = {( e RN4!/d(2,0) < 7}. (50) 


Osservazione. La misura di Lebesgue u in RN+! è invariante per le traslazioni del gruppo 
(RN+, 0), in quanto detE(t) = e'("2) = 1, ovvero la misura di Lebesgue rappresenta la 
misura di Haar del gruppo (RN+',0). Inoltre risulta 


u(B(z,7)) = u(B(0,7)) = r9*?4(B(0,1)), 
per ogni z € RN+#! e per ogni r > 0 dove 
Q+2=po+3p1+5p2+-:-+(2r+1)p+2. 


Chiameremo dimensione omogenea di RN+! il numero intero Q +2. 
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5 Costruzione del pavage intrinseco associato al gruppo 


Adesso, come in Biroli-Mosco-Tchou [28] e [29], costruiamo il pavage periodico associato al 
gruppo per il quale l’ operatore Lo è invariante. Denotiamo con B = (-1,1)N x (0,1) il 
paralleletopo di RN+!. Sia k; € Z per ognij=1,....N+1ek= (2k1,2k2,...,2kw,kw41). 
Poniamo allora 
{B.} = {ko B} (51) 

ossia il paralleletopo B traslato a sinistra di k secondo la definizione (42). 

Naturalmente non è ovvio a priori che {Bx} definito in tal modo rappresenti un pavage 
di RN+!. Comunque, con semplici argomenti, si dimostra facilmente che {Bx} è un pavage 
di RN+! e la verifica diretta di tale affermazione è resa molto semplice se si esplicitano in 
maniera più evidente le caratteristiche delle traslazioni a sinistra. «+. 

Consideriamo adesso questo pavage dilatato di un piccolo parametro € > 0, secondo le 
dilatazioni definite in (46): 


{Bi} = {D(e)Bx} = {(D(e)k)o (D(e)B)}. (52) 


Anche in questo caso è facile osservare che {Bi} rappresenta una famiglia di pavage proprio 
grazie alla (48) la quale esprime la proprietà di commutazione tra dilatazioni e traslazioni. 
Sia ora vo una funzione definita in B; possiamo definire la funzione ù come l’ estensione 
periodica di v) rispetto al pavage intrinseco {Bx}; infatti, come verificheremo tra non molto, 
posto 2' € RN+! allora esiste un unico k € (2Z)N x Z ed un unico 2 € B tali che 


. 2=koz; 
allora grazie alla proprietà delle traslazioni a sinistra 
z=k002' 


per cui definiamo 
v(2°) = do(2) = do(k7 02). | (53) 


La funzione v(-) sopra definita risulta pertanto periodica rispetto al pavage intrinseco {B,}. 
Le funzioni per le quali è soddisfatta la (53) saranno denominate funzioni Bx-periodiche. Allo 
stesso modo, fissato € > 0, se vo è una funzione definita in B*, definiamo la funzione 05 come 
l’ estensione periodica di vo rispetto al pavage intrinseco {Bi}; ancora una volta, infatti, se 
2' € RN+ allora esiste un unico k € (2Z)N x Z ed un unico 2 € B® tali che 


2' = (D(e)(k))oz 


e quindi essendo L 
z=(D(e)(k))02'= D(e)(£k71) 02, 


definiamo come prima 
v°(2') = d6(2) = d(D(e)(k7!) 0 2°). (54) 


In sostanza, come nel caso classico, la funzione v Bx-periodica definita dalla (53), quando è 
calcolata nelle variabili riscalate ponendo 


v°(2) = v(D(e7)2), (55) 


risulta essere Bf -periodica. : 
In base alla definizione di funzione Bx-periodica è immediato verificare che risulta 


Lov(:) = Lo(do(k"' 0 (-))) = (Lodo)(k7 0 (-)) (56) 


100 


ed allo stesso modo, per una funzione Bj-periodica 


Lo(v*(:)) = Lo(8(D(eV(E)"! 0 ())) = (Lo85)(DO(k1) 0 (1). (57) 


Adesso si osservi che se poniamo 


To = (1,...,tp) vettore di dim. 1 X Po 
Ti = (Cp0411---1Tpotp:) vettore di dim. 1x pi, (58) 
Tr # (Tpotpi+--+pr_1+1:-«-1£N) vettore di dim. 1 X p,, 
e successivamente indichiamo con 
Dr, = (0z11-:4%%%0) 
Da, = (02,041 si -1Izr04m ) (59) 
Dz, = (etto! se -192w) 


i gradienti nelle variabiliTo, ...,7,, allora possiamo riscrivere 1’ operatore nella forma seguente: 


q 
DO dz;(0;;(2)9:;) = divz, (a00(2) Dr, (60) 
sig=1 
dove divz, indica |’ operatore di divergenza in R° in quanto qg = po e con chiaro significato 
della notazione per aoo; inoltre si ha 


Dx, 
N Di, 
(x,BD)= YO z;b;;9.;= (0,t0B1,...,T,-1B,) : 7 
s,j=1 \ 
Di. 
= (70B;)Dz, + (£1B3)Dz, +-+ (E,-1B,)Da,; (61) 


dalle precedenti si ricava che 


L= ;» 2,,(a;;(2)d:,) + (2, BD)- 9, = 


i.j=1 
= divz,(do0(2) Ds) + (z0B1)D3, + (71B2) Da, +...t (7,-1B.)Dz, = 0. (62) 


Con questa nuova notazione risulta interessante evidenziare meglio alcune caratteristiche 
delle traslazioni a sinistra rispetto a cui Lo è invariante: consideriamo pertanto k € (2Z)NxZ, 
k= (2k1,-..:2kn,kn41) e ricordiamo che il traslato a sinistra di 2 € RN+! secondo il vettore 
k è 5 

(En,kn41) 0 (2,t)= (+ E(t)EKy;t+kw41), (63) 
dove con ky abbiamo indicato le prime N componenti del vettore & € (2Z)N x Z, mentre 


FI è semplicemente il trasposto di Kw; essendo B, per definizione, una matrice nihilpotente 
risulta che 


E(t)= exp(-tB)= >» CO gr) = 


t tr 
= Iy-tBT +3(B")+---+(-1)"3(87)". (64) 
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Con un po’ di algebra noiosa necessaria per calcolare le potenze successive di BT risulta 
infine: 


1°°%1 


DA = a ne Pa Ne ti - 
T+E(t)kv ={Z0+ko,71+F1— t(FoB;),7.+F2- t(k1B2) + 3 Fo(B1 Ba) 


= = E 
secaZrt k, “i t(k.-1B,) ca 3fr-2(B,-1B,) as 
SY 


r! 


+ 


t"Ko(B1(B2(:--(B.-1B.)--))))}, (65) 


dove 
ky = (Kos F1,...,k7) = (2k1, 2ka,...,2kw) 


in consistenza con la notazione introdotta tramite le (58). 


6 Definizione del problema 


Sia G un dominio limitato di RN#!; in particolare supporremo G di tipo cilindrico, ossia 
G=Qx[0, T) dove Q è un aperto limitato di RN. In G consideriamo tin particolare problema 
di Dirichlet ai valori al contorno associato ad un operatore Li, (di cui tra poco chiariremo 
il significato) ed i cui coefficienti sono rapidamente oscillanti, o meglio sono Bi.-periodici 
secondo la definizione (3.19). i 

Indichiamo con Sx=p l’ insieme di tutte le funzioni a : RNH#! —; M?*, con a(2) = 


(a;;(2)) = (4;;(2)) per ogni 2 € R#! e per ogni i,j= 1,.. -,9, tali che: 
A1) ai; è Bx-periodica per ogni i,j=1,.. RI 

Az) a;; € L°(RNH); 

A3) esiste una costante 4 > 0 tale che 


9 9 9 
29% < &;(2)&6 < pi CE 


Hizi ij=1 i=1 
per quasi ogni 2 € R+! e per ogni € € R9. 


Alle funzioni a che godono delle proprietà appena citate, associamo la seguente famiglia di 
operatori 


L= = ds; (0°(2)d;,) + (2, BD)-d (66) 
dove abbiamo posto 
a°(2) = a(D(e71)2) (67) 


con £ piccolo parametro positivo. Si osservi, in consistenza con la (55), che le funzioni di; 
risultano Bf-periodiche in BN+!. A questo punto un approccio diretto'alla formulazione debole 
del problema di Dirichlet associato all’ operatore indicato nella (66) risulta immediatamente 
insoddisfacente, se condotto con i metodi standard della teoria degli operatori parabolici 
astratti largamente discussa da Lions. 

Allora un modo per superare questa difficoltà, è quello di approssimare il problema origi- 
nale con una sequenza di problemi variazionali ben posti che hanno soluzioni uniche. L’ idea 
fondamentale di questo lavoro è stata quella di effettuare un particolare procedimento di rego- 
larizzazione per la classe di operatori ultraparabolici studiata: l’ approssimazione effettuata 
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I abbiamo definita regolarizzazione subellittica. Definiamo allora il regolarizzato subellittico 
dell’ operatore L° tramite la seguente 


suo = eV? + L* (68) 


dove la parte regolarizzante è costituita dal termine e°Y?. Si osservi che la potenza con cui 
compare il parametro e non è casuale ma strettamente legata, come vedremo, al procedimento 
di omogeneizzazione che applicheremo a questi operatori. In sostanza questa regolarizzazione 
ci consente di approssimare l’ operatore originale, che è un operatore ultraparabolico, con 
una sequenza di operatori che godono di proprietà molto potenti ed interessanti. 

Infatti la classe di operatori del tipo (68), per e > 0 fissato, verifica stime subellittiche 
(si veda [27]), da cui il nome scelto per il procedimento di regolarizzazione; inoltre grazie 
alla (H;), sempre per e > 0 fissato, l’ operatore Lî,; risulta coercivo, come verificheremo 
tra non molto, in virtù della disuguaglianza di Poincaré valida per campi vettoriali di tipo 
Hòrmander (per questo risultato si rimanda al lavoro di Jerison [30]}). Pertanto la classe 
di operatori definita nella (68) ci consente in modo agevole di lavorare nel contesto delle 
soluzioni deboli sfruttando, come vedremo, il lemma di Lax-Milgram grazie alla sopraccitata 
validità della disuguaglianza di Poincaré per campi vettoriali di tipo Hòrmander. 

Infine, sottolineiamo che la regolarizzazione effettuata è tale da preservare 1’ invarianza 
per traslazioni a sinistra definite nella (42) come enunciato nella seguente proposizione di 
verifica immediata 


Proposizione 7 Ogni operatore 
Lo,sub = ey? + Lo, 


dove Lo rappresenta l’ operatore a coefficienti costanti definito nella (40), è invariante per le 
traslazioni a sinistra di (RN+ 0), ovvero 


Lo,sub(u(2' 0z))= (Lo,sub)(2" o 2) 
per ogni z € RN+! con 2' € RN+! fissato. 
Adesso siamo pronti alla formulazione esatta del problema che intendiamo affrontare: per 
un fissato e > 0, consideriamo il seguente problema di Dirichlet ai valori al contorno su G 
{ eY?u+L"u = finG, 


u“ lac = 0su 0G, (69) 


dove f è, al solito, una funzione assegnata su G; in particolare, nell’ ambito della teoria dell’ 
omogeneizzazione, ci interessa stabilire il comportamento di u° per € che tende a zero. Prima 
di passare alla formulazione debole del problema (69) abbiamo bisogno di specificare in quali 
spazi funzionali andiamo a cercare le nostre soluzioni deboli. Allora denotiamo con H#(G) lo 
spazio di Sobolev degenere associato ai campi vettoriali indicati nella (H}) cioè 


H(G)={u € L°(G)/u., € L°(G)Vj=1,...,g,YuE L?(G)}, — (70) 
fornito della norma naturale 
lu llg=ll Ilio + Il Das Iwa + Yu [SE (71) 
facciamo presente che in questi spazi vale un teorema tipo Meyers-Serrin, cioè 
C°(G)N H(G) è denso in H(G) 


(si vedano [31], [32], [33]). Pertanto è naturale denotare con Ho(G) la chiusura di C5°(G) in 

H(G). Di fondamentale importanza risulta poi il seguente teorema che discende da quello più 

generale provato in [30] riguardante, appunto, campi vettoriali che soddisfano la condizione 
‘ di Hòrmander 
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Teorema 8 Esiste una costante geometrica c > 0 tale che 


i |u|?dz < c( / |Dz,u|?dz + / IVu]?dz) (72) 
o o o 
per ogni w € Ho(©) e per ogni aperto limitato O di RN+. 


In sostanza, il teorema appena ricordato asserisce la validità della disuguaglianza di 
Poincaré per i nostri campi vettoriali. In particolare, segue che la norma naturale definita 
nella (71) e la seguente norma 


lu |ioo=Il Dart lio + | Yu Il) (73) 


sono norme equivalenti nello spazio Ho(0). Per cui d’ ora in avanti, salvo diversa indicazione, 
quando ci riferiremo allo spazio #o(0) lo supporremo sempre fornito della norma || - |lxx(0) 
definita nella (73). Possiamo ora formulare il pricipale risultato riguardante il problema di 
Dirichlet per l’ operatore L5,, in G. 


Teorema 9 Sia f € L?(G). Allora esiste una ed una sola u° € Ho(G) che risolve il problema 
di Dirichlet per e > 0 fissato 
(P.) { Li up0 f inG, 


ulae = 0su0G 


nel senso che 


{ M°(u,v)=-fsf(2)v(2)d2 per ogni v € Ho(G), (74) 


u E Ho(G) 
dove M° : Ho(G) x Ho(G) + R è la seguente forma bilineare 


M°(u, v)= f(@Daw, Ds,w)dz+e? | YuYodz+ / u(Yv)dz (75) 
Te (e G 


La verifica è immediata tramite I’ applicazione del lemma di Lax-Milgram. 
Riconsideriamo il problema (69) nella seguente forma 


Lit = f inG, 
{ u € Ho(6) VAI 


dove supponiamo f € L?(G). Consideriamo adesso la successione (u“) di soluzioni del prob- 
lema (76) al variare del parametro e. È facile verificare che 


1 \ 
E || Yu © ta | Da [PERSE S Ilex)l] « |) (77) 


da cui 
ce luo f ll] 4 lx), (78) 


dove c, è una costante che dipende dal parametro e. Il passo successivo è stato quello di cercare 
stime a priori su u°, indipendenti dal parametro e, in uno spazio funzionale più grande dello 
spazio Ho(G) e sicuramente meno regolare di questo, ma comunque più congeniale ai nostri 
scopi. Definiamo allora 


Hî°(G) = {ue L?(G)/ us, € L’(G)Vj=1,...,9} (79) 


fornito della norma 
lula=il Pat Io + Il è II (80) 
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e si osservi che 
C°(G)NHF°(G) è denso in H7°(G) 
sempre per la validità di un teorema tipo Meyers-Serrin in questi spazi (vedi [31], [32], [33]); 
quindi è naturale definire con H#7°(G) la chiusura di C5°(G) in H"°(G). 
Con argomenti standard di tipo classico si arriva facilmente a verificare la validità. delle 
seguenti 


Il u° PETE Cri (81) 
| u° [IL (@)£ c2 (82) 

ed 
€ || Yu® |l22(6)£ ca (83) 


dove c1, ca, c3 sono delle costanti reali positive che non dipendono dal parametro e. Grazie alla 
(81) ricaviamo che è possibile estrarre una sottosuccessione, ancora denominata allo stesso 


modo, tale che per e + 0 d. 
u + u in H5°(G) - debole . (84) 


7. Costruzione dell’ operatore omogeneizzato 


Introduciamo in questo paragrafo il metodo dello sviluppo asintotico di cui abbiamo già 
discusso in precedenza. Al solito, cerchiamo uno sviluppo asintotico per u“(z), come funzione 
del parametro e: ipotizziamo in modo del tutto euristico che u" ammetta uno sviluppo in 
due scale della forma 


u°(2) = vo(z, D(e71)2) + e'u1 (2, D(e7!)2) + e2u2(2, D(e7")2) + a (85) 


dove le funzioni w;(z,() sono funzioni Bx-periodiche nella variabile C (che chiameremo vari- 
abile microscopica) per ogni z € G (che chiameremo variabile macroscopica). Quindi noi 
postuliamo l’ esistenza di funzioni regolari w;(2,6) definite in G x RN+!, indipendenti da €, 
B,-periodiche nella variabile € e tali che per € = (D(e7')2) il lato destro della (85) sia uno 
sviluppo asintotico di u"(2). 

Facciamo notare che quando consideriamo le funzioni 


u;(2, (D(e7!)2)) (86) 


queste risultano essere invece Bj -periodiche nella variabile microscopica. Esattamente come 
nel caso classico, si procede inserendo la (85) nella (69) ed identificando le potenze successive 
di €; si noti che è proprio in questa identificazione di potenze che gioca un ruolo fondamentale 
I’ esponente con cui appare il parametro € nella parte regolarizante della (68). Anzitutto 
ricordiamo che in consistenza con la notazione (58) possiamo scrivere 


(D(£)712) = (e7170,e78#1,6572,. e VE, 614) = 


ni (€01 E En) sE siC T). (87) 
Pertanto quando poniamo ( = (D(£7')2) intendiamo, în base alla (87), 
= (Enia cafatl= (€ 7). (88) 


Assegnata poi una funzione regolare ®(2,6) di due variabili, definiamo la funzione 


®*(2) = D(2, D(e !)z) (89) 
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ed osserviamo che 1 
Dz,9"(2) = (D=,® + Pr9(2, D(e-1)2) 
Da,®*(2) = (Dz,® + 5Dx,9)(£, D(e1)A) 
Dz,9°(2)=(Dz,®+ 30,9, D(e71)z) 
1 
Dz,®°(2) = (Dz,D+ pari Di, ®)(2, D(e!)2) 


4,8°(2) = (,0+ 59,8)(2, D(e-1)2). 


Facciamo notare che in questo contesto stiamo considerando le variabili z e C indipendenti, 
per cui i corrispondenti operatori di derivazione commutano tra di loro. Adesso, con questo 
in mente, analizzando il singolo comportamento dei tre termini che compongono l’ operatore 
Liub quando sono applicati a funzioni del tipo (89), otteniamo che: 


Liub®" = [(E72(Lo + Y?,) +67 (£1) +6%(£2 + 2YHe,7))®](2, D(e71)2) (90) 
dove per G = (z, D(£-1)z) gli operatori sopra indicati sono i seguenti 
q 
Lo= )7 de;(0;;(6)%,) + (€, BDe) — d, (91) 
i.j=1 

q q 

Li = )7 de, (0;;(0)%,) + DO d(0:;(0)9:,) (92) 
î.j=i i,j=1 

L2= IO ds,(0;(0)0:,)+ (7, BD.)- è. (93) 


t.g=1 
Pertanto utilizzando la (85) e la (90), I’ equazione (69) è soddisfatta, sotto l’ ipotesi che la 
matrice a e le funzioni u;(z,C) siano sufficientemente regolari, se 


(Loto +YE, uo)(2, D(e7!)z)=0 inG, (94) 
(Lou, + Yu + Lito)(2,D(e7!)2)=0 inG, (95) 
(Loua + Ye, u2 + Li, + Lato + 2Y Yes tto)(2, D(e7!)z) = f inG. (96) 


Mostriamo adesso come dalle (94), (95) e (96) si possa costruire 1’ operatore omogeneiz- 
zato. Procedendo come nel caso classico, determiniamo le funzioni u; risolvendo i problemi 
seguenti nei quali consideriamo z come un parametro fissato e nei quali le condizioni di 
By-periodicità svolgono il ruolo di condizioni al contorno: 


1 (Covo + Ye, to)(z,-)=0 in B, 07) 
Uo(z,-) Bx-periodica ; 
(Low +Y#,)(z,:)=—(L1u0)(z,-:) in B, (98) 
u;(2,-) Bx-periodica ; 
(Low + Y?_u2)(2,-) = f(2) — (£200 + 2YerYzeto+ L10;)(2,-)in B, (99) 
u2(2,-) Bx-periodica . 
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Consideriamo questi problemi nel contesto delle soluzioni deboli. Denotiamo allora con 


L»,,(B) la chiusura dello spazio C$,(RN+!) rispetto alla norma || - |lz=xa), dove 
.Ce.(RN+) = {ue C°(RN+!), u By-periodica}. (100) 


Allo stesso modo definiamo #,»-(B) la chiusura dello spazio C®,(R+!) rispetto alla norma 


per 
Il - Ila) che è una norma su C$2.(R+!) proprio a causa della periodicità. Inoltre indichiamo 


con Hy(B) il sottoinsieme di Hper(B) di tutte le funzioni con valore medio nullo su B, ovvero 


Hg (B) = {u € Hoer(B)/7g7 f, (Od =0). (101) 


A questo punto possiamo introdurre un risultato molto importante di esistenza di soluzioni 
per un problema ai valori al contorno sulla cella di base B associato ad operatori subellittici. 
Lemma 10 Sia F € L},,(B). Allora il problema 


Ad, Ù) su e fa Fd per ogni 0) € Hper(B) , (102) 
$ E Hy(B), 


dove 
A($, 4) : Hper(B) X Hper(B) + R 


è la forma bilineare definita da 
q 
A@)=L fr ass ODE ADa dk + fVerdte dc + / WMA (108) 


ammette una unica soluzione se e solo se 


L F(()dC=0. (104) 
Vorrei solo sottolineare a tal proposito l’ importanza delle seguenti relazioni 
f,Aen(Vdc — 0 per ogni j = 1,...,9 e per ogni € Hper(B) , (105) 
e 
/ Vu(0)d(=0 per ogni u € Hper(B) (106) 


che non solo ci consentono di passare alla formulazione debole dei problemi (97), (98) e (99) 
grazie alla validità di alcune formule di integrazioni per parti ma, come vedremo, risultano 
fodamentali anche per varie considerazioni relative allo studio dei problemi appena citati sulla 
cella fondamentale. ° 
Consideriamo adesso il seguente lemma, caratterizzante le proprietà della estensione periodica 
rispetto al pavage intrinseco {Bx} della soluzione del problema (102): 


Lemma 11 Per F € L?(B), a media nulla su B, si consideri l’ unica soluzione debole $ € 
Hy(B) del problema (102); allora l’ estensione periodica $ di $ rispetto al pavage intrinseco 
{Bx} è una soluzione, nel senso delle distribuzioni, del problema 


Lid+Y},6=f ing, (107) 


. dove F rappresenta l’ estensione periodica, rispetto al pavage {Bx}, di F. 
pp 9g 
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Concludiamo, per finire, con un lemma caratterizzante le proprietà dell’ operatore Lo+Y}, 
sotto ipotesi di maggiore regolarità sui coefficienti. 


Lemma 12 Supponiamo che l’ ipotesi (A) sui coefficienti d;j(2) € Sxrp sia sostituita da 
quella più forte a;;(2) € C°(RN+!), ferme restando le altre ipotesi. Allora l’ operatore 


Lo+Y}, i (108) 
è ipoellittico. 


Dimostrazione. L’ operatore è ipoellittico come si può vedere utilizzando le tecniche di Kohn- 
Radkevic (vedi [34], [35]). Per una formulazione esplicita del risultato si veda M. E. Taylor 
[36], cap. 15, Teorema 1.2. 

Adesso, procedendo formalmente come nel caso classico, grazie al Lemma 10 possiamo 
determinare le soluzioni dei problemi (97), (98) e (99). 


Studio del problema (97). 
Il problema (97) può essere formulato in modo debole come segue: 


{ A(v0,9)=0 Ve Hper(B), (109) 


Uo(2,:) E Hper(B) 


dove A è la forma bilineare introdotta nella (103). Ricordiamo che in questo contesto 2 gioca 
il ruolo di un parametro. Allora grazie al Lemma 10 concludiamo che uo(2, +) è determinata 
dall’ identità integrale presente nel problema (109) a meno di una costante. Dalla stessa 
identità, prendendo % = uo(z, +) e sfruttando la coercività della forma bilineare A, si ricava 
che 

1 Il vo(2,*) Il,.-(8)S A(t0, vo) = 0 (110) 


ovvero che uo(2, -)=costante, quindi 


uo(2,:) = uo(2) (111) 


Studio del problema (98). 
Grazie alla (111) la formulazione debole del problema (98) diventa: trovare u1(z,+) € Hper(B) 
tale che 


: duo 
A(,4)=- D dee fas Deva VW € Hpe-(B). (112) 
Allora consideriamo per j = 1,...,g il problema 
A(wi, 4) = — fg(a(6)e;, Dr W(0))dC VW E Hy(B), (113) 
wi € Hy(B). 


Osservato che l’ applicazione Fi : Hy(B) + R definita tramite la seguente 
va F()=- f (06; DO) 


è un funzionale lineare e continuo su Hy(B) e ricordato che in virtù della relazione (105) 


risulta x 5 
Qi \ de = 
= lea =o, 
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sempre dal Lemma 10 discende che esiste ed è unica la soluzione del problema (113) per ogni 
j=1,...,9. Quindi una soluzione generale del problema (112) è 


Dale, 


LOGO ci(2) (114) 


dove c1(2) è una costante additiva dipendente, quindi, dal. parametro 2. 

Studio del problema (99). 

Nel problema (99) chiaramente |’ incognita è rappresentata dalla funzione w2 essendo, al 
solito, z un parametro. Pertanto grazie al Lemma 10, u2 esiste ed è unica se e solo se 


f,(Cao + Cdk = f, 1(A)d< = IBINA), (115) 
dove |B| indica la misura di B; si noti che nella precedente abbiamo sfruttato il fatto che 


YaeYe,rUo(2) = YerYasto(2) =0 


in quanto i due campi Y;x e Yer, rispettivamente nelle variabili macroscopiche e micro- 
scopiche, commutano tra di loro. La condizione (115) rappresenta, in analogia col caso 
classico, proprio !’ equazione omogeneizzata che stiamo cercando: infatti osservato che 


CI 


Liu, = LE P (as(6 + x TAGLIO OT), 


in 


si ha, utilizzando la (114) e la proprietà (105) 


do, 
fr (fade = © Fal Lf, Dau(03 sta 
+ x ha EA: (Dan MITE = > ha DE IO an (0)d)ZI = 
a? Uo 
"2 È Dan) 7E De (Od, (116) 
per cui la (115) si esprime come 
3 d,.(ate"9,,u)+Yu=f inG, (117) 
dove per î,j = 1,...,9 
om _ a 
ai; = NJ, (a;;(6) )+ Do 0u(0 )dc. (118) 


L’ equazione (117) è l’ equazione vmopenzizzala mentre i coefficienti af?” sono i coeffici- 
enti omogeneizzati. A questo punto un po’ di considerazioni sono necessarie. In sostanza 
il procedimento di regolarizzazione subellittica presenta un duplice vantaggio: anzitutto, 
come abbiamo visto in precedenza, ci consente di lavorare nel contesto delle soluzioni deboli 
sfruttando a tal proposito il Teorema 9 e per i problemi ai valori al contorno sulla cella di 
base B associati ad operatori subellittici il Lemma 10; inoltre, ipotizzando uno sviluppo della 
forma (85), ci consente (per il momento formalmente) di costruire I’ equazione omogeneizzata 
‘ (117), la cui forma è, però, quella di una equazione di Kolmogorov-Fokker-Planck. Infine è 
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immediato verificare che, nelle nostre ipotesi, la matrice dei coefficienti omogeneizzati risulta 
simmetrica. 

Nel prossimo paragrafo mostreremo un particolare metodo di regolarizzazione per funzioni 
u € L3.,(B) che preserva la periodicità rispetto al pavage intrinseco {Bx}. Sfruttando questo 
metodo per approssimare i coefficienti degli operatori di interesse con coefficienti. regolari, 
vedremo che sarà possibile costruire una classe di funzioni test di fondamentale importanza 
per la prova del risultato di convergenza. 


classe di funzioni test 


Intendiamo ora approssimare i coefficienti del nostro operatore con coefficienti regolari che 
rispettino la periodicità rispetto al pavage intrinseco. È naturale quindi utilizzare convoluzioni 
che rispettino la struttura di gruppo con dilatazioni. Un simile approccio al problema della 
regolarizzazione dei coefficienti, per equazioni tipo Kolmogorov-Fokker-Planck, appare anche 
in [25], [16], [26], in forma leggermente diversa, essendo formulato in una diversa prospettiva. 
Pertanto, sia w una funzione non negativa, a valori reali appartenente a C5°(RN+!) ed avente 
le seguenti proprietà: 

1) w(2) = 0 se |||z]]| > 1, dove con |||: ||| ci riferiamo alla norma introdotta nella Definizione 

, 
2) Sana W(2)dz = 1. 
Adesso sia 7 > 0: allora la funzione 


cn (2) = é(D(n°!)2)p(9+2 (119) 


è non negativa, appartiene a C$° (RN+1) e soddisfa le seguenti: 
a) wy()=0se [Ilz]l>n,e 
guy Wp (2)dz = 1; 
la proprietà (a) discende, in particolare, dalla (N) della Proposizione 5. Sia ora u(z) una fun- 
zione By-periodica; supponiamo inoltre che u(2) € Lî,c(R+!) e definiamo K-regolarizzazione 


di u(z) la seguente funzione 


un(2)= fi. MA) 2A? ; (120) 


Osserviamo che, nelle nostre ipotesi, tale definizione è ben posta; infatti essendo w,(z) a 
supporto compatto si ha che 


u(2')wn((2)7! 0 2)|d2' < suplw,l |u(2')|d2" < 
(IA, bici 7 Jnen-sesii<o 


<p (042) lu(2)|d2" (121) 
III(=")-toz]II<m 


e quindi grazie alla disuguaglianza di Hélder risulta 


- (042) ul2)ld2 < n-(242) 2,9) u(2')|?d2)} = 
1949 prc EN STUB MA lle) 
= 28 u(B(0,1)}A(f  pule)tez)). (122) 
B(2,n) 


dove con B(2, 7) intendiamo l’ insieme indicato nella (50). 
Le proprietà della funzione u,(z) sono riassunte dal seguente 
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Teorema 13 Sia u, la funzione definita dalla (120) con u(2) ew,(z) che soddisfano le ipotesi 
precedentemente fatte; allora sono verificate le seguenti: 

M;) u,(z) è By-periodica, 

Ma) un(2) € C©(RN#), 

Me) un(2) € L2,,(B), || ur(2) IlnuS<Il (2) Îlzuo € 

limp0 || un(2) — (2) ||z=x)= 0 per ogni compatto I € RN+!. 


La (M;) discende facilmente dalla definizione di By-periodicità, mentre per la (M:) ed (M3) 
si procede come nella teoria dei mollificatori classici (vedi, ad esempio, Adams [26], Teorema 
2:18); 

Adesso consideriamo il X-regolarizzato, secondo la definizione (120), dei coefficienti a;;(2) 


a" (2) = (al(2) = (f.,, 05 (#0, ((#)-! 0 2)d2) (123) 


dove a;;(2), per i, j = 1, ...,g, soddisfano le ipotesi (A1), (A2) ed (43); 

in particolare essendo a;;(2) € L°(RN+!), deduciamo facilmente che a;;(z) € L7,.(RN+'). 
Pertanto dal Teorema 13 ricaviamo che i coefficienti X-regolarizzati soddisfano le seguenti 
proprietà: 

Ai) ai;(2) è Bx-periodica per ogni i,j = 1,...,9; 

A5) aîi(2) € C®(R"H), at.(2) € LE,,(B); limn=0 || 08;(2) — a4;(2) Il: per ogni compatto 
K di RN+, quindi limo || a7;(2) — 4;;(2) |lza,.(ay= 0 per ogni i,j= 1,...,9; 


‘per 


A}) esiste una costante y > 0 (la stessa introdotta nella (A3)) tale che 


q 


g q 
a < D ag;(2)&6; < ADE 


i.j=1 
per ogni 2 € RN+! e per (£1,...,£) € R® (si osservi che in questa situazione le variabili £; 
non rappresentano variabili microscopiche periodiche ma le componenti di un generico punto 


di R°). 
Naturalmente lo stesso procedimento di regolarizzazione può essere applicato, con una lieve 
modifica, anche a funzioni Bi-periodiche; posto infatti per e > 0 fissato 


a; = (fi) 0 2a}; (2)d2) = 


= (Li wy((D(e7)21)7! o D(£7!)2)a;;(D(e)2/)dz') (124) 


è immediato verificare la validità di proprietà, analoghe a quelle appena enumerate per i 
coefficienti a?; (2) tenendo però presente che adesso i coefficienti della matrice a; (2) risultano 
periodici rispetto al pavage intrinseco {Bf}. i 

Analizziamo adesso il comportamento della soluzione del problema (P°), pere > 0 fissato, 
con coefficienti X-regolarizzati secondo lo schema appena introdotto; quindi consideriamo il 
seguente problema 


Liwu,=f inG, 
sui 12 
{ ui € Ho(6), VERI 
in cui intendiamo a 
sm=eY?+ DO 0e;({a5(2));;0,,) +Y (126) 
i.j=1 


dove la matrice dei coefficienti è la X-regolarizzata della matrice a"(2) secondo il procedi- 
‘ mento indicato nella (124). Si osservi che in tutte le considerazioni che seguiranno in questo 
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paragrafo consideriamo il parametro e > 0 fissato. Grazie alla proprietà (43) è immediato 
verificare che esiste una ed una sola ur € Ho(G) che risolve il problema (125) per n>0 
fissato, nel senso che 


Mi(u0)=-faf(2)v(2)dz Vu e Ho(G), 
{ ut € Ho(0), NS 


dove Mj : Ho(G) x Ho(G) + R è la seguente forma bilineare 
Mi (u,v) = L (at Dz,u, Dz,v)dz + e? È YuYudz + b u(Yv)dz. (128) 


Il comportamento della successione di soluzioni (u5) al variare del parametro 7 è stabilito 
dalla seguente 


Proposizione 14 Siano u5 € Ho(G) ed u € Ho(G) le uniche soluzioni deboli rispettiva- 
mente dei problemi (125) e (76); allora, pern +0 


un + u in Ho(G)-debole; (129) 
in particolare essendo Ho(G) + L?(G) in modo compatto segue anche 
u, +u" in L?(G) -forte. Po (130) 


Inoltre è importante far notare che sono facilmente ottenibili stime del tipo (81) e (82) per 
la soluzione «5 del problema (125), ovvero 


| U Wly,zo G < Ci (131) 
0° (6) 


Il us Ilz26)< c2 (132) 


dove le costanti c1 e c, sono indipendenti sia dal parametro n che dal parametro e. 
Premesso ciò ritorniamo al nostro obiettivo fondamentale che è quello di costruire una 
particolare classe di funzioni test. Poniamo allora 


vr (2) = vo(2) +e81 (2,6) + e282(2,0) (133) 


per € = (D(e7!)2) dove vo(2) E C$°(G) mentre 07(z,-), per j = 1,2, rappresentano le 
estensioni periodiche rispetto al pavage intrinseco {Bx} delle funzioni v7(2,-) che adesso 
specificheremo: 1 

Ci) vi(2,0) = YI, Bs(C)PEla, nella quale 95(6) per s = 1,...,4 è l’ unica soluzione, in 
senso debole, del problema (notare che d’ ora in avanti, per non appesantire la notazione, 


sottointendiamo il simbolo di sommatoria sugli indici ripetuti) 


(Pî) ELI Ù > * (134) 


dove 
A" :Hy(B) x Ha(B) HR 


è la forma bilineare definita da 


AGI = [LOUIE f Vesdk vd - [ esa. (135) 
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L’ esistenza e l’ unicità di una soluzione debole del problema (P7), per 7 > 0 fissato, grazie 
alla (45), è garantita dal Lemma 10; inoltre con discorsi analoghi a quelli utilizzati per provare 
la Proposizione 14, si verifica immediatamente che per 7+ 0, 


Bi +° in Hy(B)-debole (136) 
per ogni s= 1,...,g dove 8*(C), per s= 1,...,9, è I’ unica soluzione debole del problema 
A(p°,%) = - Ja Gis (0) de; W(C)de, (137) 
B° € Ha(B), 


in cui A è la forma bilineare definita sempre dalla (135), considerando in questa i coefficienti 
a;;(6) anzicché quelli K-regolarizzati. Facciamo osservare che il problema (Pf) rappresenta 


la formulazione debole del problema aggiunto sulla cella di base B, per s= 1,...,9 
(LI) B:() +VE A; (0) = - Edel (0) in B, js 
B;(-) B, - periodica. 


Infatti, grazie alla simiarita dei coefficienti, l’ aggiunto dell’ operatore con coefficienti X- 
regolarizzati (£ +Y},) risulta semplicemente l’ operatore ((£3)* +Y?,) dove 


LÀ 
(£8)° = DO de.(ai;(0)A,) — (€, BDe) + dr. 
ij=1 
Grazie alla regolarità dei coefficienti e dei dati, grazie al Lemma 11 ed al Lemma 12 (ipoel- 
liticità dell’ operatore considerato) si verifica facilmente che 1’ estensione periodica rispetto 
alla variabile microscopica di periodicità 


LI dvo(z 
=Lh (eeo. (139) 
è una soluzione di 
((£L3)*01 + Yi, 01 +(£L7)*vo)(2,:)=0 in RN (140) 


dove (£7)* = (£1), ovvero l’ operatore £1, i cui coefficienti sono X-regolarizzati, è simmetrico. 
Ca) v3(2,6) = Din (MA dove x!(6) per 1,s = 1,...,9 è l’ unica soluzione del 


seguente problema aggiunto BS cella elementare 
mv) = - Sala" bi, (CdC Vy € Ha(B), 
(Pi) {4 x cHy (B), - n (141) 
in cui A” è la forma bilineare definita dalla (135) mentre per ogni /,5= 1,...,9 
alert = TE (O + (0A (142) 
e 
br, (6) = an, (6) + at; (09:95 () + de; (ai (6)Ar(6)) (143) 


Si osservi che i coefficienti indicati nella (142) rappresentano i coefficienti dell’ operatore 
omogeneizzato 


(KP), = i» 9: (aff” "9;) — (2, BD) ca è, (144) 


tj=1 
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relativo all’ operatore (Lib) e che possono essere facilmente costruiti seguendo il procedi- 
mento formale indicato nel paragrafo (7). Anche in tal caso si può verificare che la simmetria 
di (a7;(C)) implica quella di (apra), 

D’ altronde costruendo i coefficienti dell’ operatore (Là°M), relativo però allo operatore (L57), 


avremmo ottenuto per questi, per (;a=1,.. -39 
1 
prom — mJ ar, (6) + af; de, wi (C))d 145 
I |B| R ts (6) Dei a (0) 6 ( ) 
nella quale w;(6) per s= 1,...,gè l’ unica soluzione del problema 


{ B" (wp, 6) = — Sp al.(VAP(C)d Vé € Hg(B), (146) 
Wi (SÌ Hy(B), 


dove 
B" : H4(B)x Hy(B) 4 R 


è la forma bilineare definita da 


B"(0,0)= [VA AMA [Verte 4A + forervic. (147) 
Ma, come si può facilmente verificare con semplici argomenti, risulta per s = 1,...,9 
Bi(6) = wi(C) quo. in B. (148) 


Pertanto l’ operatore indicato nella (144) rappresenta |’ aggiunto dello operatore (Lf°”) 
(LE) = 04 apt PA 


ossia l’ omogeneizzato dell’ aggiunto dell’ operatore (Lil) è proprio l’ aggiunto dell’ operatore 


omogeneizzato (Lh°") grazie alla simmetria dei coefficienti algra e grazie alla validità della 
(148). 

Chiusa adesso questa parentesi necessaria per evitare eventuali ambiguità, osserviamo che 
per 7 > 0 fissato l’ esistenza e l’ unicità di una soluzione del problema (P7) sono garantite 
grazie alla proprietà (45) ed al Lemma 10: 

in particolare, per l’ applicabilità del Lemma 10, si tenga presente che per ogni /,8= 1,.. @ 


fata — in ac= o 
in quanto | 
fa (a (0s 0) =o 


grazie alla (105). 
Di fondamentale importanza risulta poi la seguente 


Proposizione 15 Indichiamo con Afer la matrice q x q i cui coefficienti, costanti, sono 
indicati nella (142) e sia poi B" la matrice qX q definita nel seguente modo 


(B?),, = TÀ bi (C)W(C)d6 per l,s=1,...,9 (149) 


dove 4(C) è una arbitraria funzione di Hy(B) ed i coefficienti b,, per 1,5=1,.. -,9, sono 
quelli indicati nella (143). Allora 


Pi — 4° pern-+0 (150) 
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B" >-B pern->0 (151) 
dove i coefficienti di A*°" sono dati da 


abem = f[ (an (6) + an (0A,8°()d, 


mentre la matrice B è 


(Bu = fu (OVIL 


nella quale, per l,s= 1,...,9 


bus (€) = 01» (6) + a (0): 8" (6) + de (an (0)P° (0). 


Osserviamo che grazie alla Proposizione 15, procedendo come pet la verifica della Propo- 


sizione 14, si ottiene che per ogni /,5= 1,...,9 
x + x" in Xy(B)-debole (152) 
dove x! è l’ unica soluzione per /,s = 1,...,q del problema 
A(x°,9) = a Se (aper (aa dis (0)(01(@ (153) 
x € Ha(B), | 


Sottolineiamo che il problema (P?) rappresenta la formulazione variazionale del seguente 
problema aggiunto sulla cella di base B, per ,5= 1,...,9 


(EA) + Yet (0) = (al? — bi, (0)) in B, 
i lu i Ù (154) 


dove i coefficienti af°”" e b7, (C) sono indicati rispettivamente nella (142) e (143). Come prima, 
grazie alla regolarità dei coefficienti e dei dati, grazie al Lemma 11 e 12, si può mostrare che 
l’ estensione periodica 


86,0 = D WOZO (155) 
è una soluzione di 
((£2)°37 + Y2,02+ (L9)"07 + (CI) "vo)(2,:) = (L8°7)vo(2) in RNH, (156) 


dove l’ aggiunto dell’ operatore con coefficienti K-regolarizzati (£1) è 
q 
“i }; dx;(a7;(6)9z,) Toi (x, BD) “È à, 
sg=i 


mentre (Là°")" rappresenta l’ operatore indicato nella (144). 
Inoltre, grazie alla proprietà di convergenza espressa dalla (150), risulta per 7-+ 0 


(LR9)"9(2) + (L*)g(2) (157) 


per ogni g(2) € C°(G), con 


(£hom)° = D a. (ata, )- (2,BD)+d (158) 


i.j=1 
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dove afom rappresentano i coefficienti della matrice A*°" o, se si vuole, i coefficienti definiti 
nella (118) come si è fatto osservare in precedenza. 

Quindi ricapitolando abbiamo definito la funzione vi (2) tramite la (133) dove, per j = 1,2 
dj (2,:) € CS,(RN+) per ogni 2 € G e 07 (-, 6) € CS°(G) per ogni ( € RN#!, Di conseguenza 


vi(2) € C®(G) (159) 


e naturalmente 


vr(2) € Ho(G). (160) 


A questo punto grazie alla regolarità dei coefficienti e della funzione v£ 
direttamente il seguente sviluppo aggiunto 


(Liu6)"v7(2) = 


= (Lio) vo(2) + (Li) 01 (2, D(e-!)2) + e°(L50)"82(2, D(e-)2); (161) 


ricordiamo, a tal proposito, che si era definita la funzione v;(2) tramite la (133) dove co- 
munque intendevamo ( = D(e7!)z. 
Allora tenendo presente la (90) è immediato verificare che per j = 1,2 


(2) possiamo calcolare 


(Liub)"07(2, D(e7!)2) = e72{((£9)° + YA, )0M](2, DE) A+ 


+ e I(C1)*37](2, D(e!)2) + e°I((L9)" + 2Y..1Ve,+)07](2, D(E-1)2). (162) 
Pertanto, grazie alla precedente, lo sviluppo aggiunto indicato nella (161) può essere espresso 
nella forma 
(Liu6)"v5 (2) = 


= e ![((L6)" + Y#,)07 + (LI) "vo](2,0)+ 
+e°I((L8)° + Y?,)82 + (L1)"37 + (2) *vol(2,0+ 
+e! {((£2)° + 2Y2,4Ye,7)0? + (£1)"#2](2,0)+ 
+e2[((£2)" + 2Yz:Y,7)1 + Y2,00](2,0)+ 
+ e°1Y201](2,0) + e4Y2,02](2,0) (163) 
dove abbiamo nuovamente esplicitato la presenza della variabile microscopica di periodicità 


Gc 

Adesso dovrebbe essere chiara la particolare scelta effettuata sulle nostre funzioni; si ricordi 
infatti che 07 (2, C) definita dalla (139) è stata costruita in modo da soddisfare la (140) mentre 
02(2, 6) definita dalla (155) soddisfa invece la (156). i 

Queste proprietà ci consentono di poter riscrivere il precedente sviluppo aggiunto nella 


seguente forma 
(Lib) "Un (2) = (£5°7)"v0(2) +697 (2,6) +e%92(2,0) +e°99(2,0) +e‘92(2,0) (164) 


dove con g}, per j = 1,...,4, abbiamo indicato i termini racchiusi entro parentesi quadre 
nella (163) che moltiplicano potenze successive di e. Nel prossimo paragrafo mostreremo come 
sfruttando questa particolare forma di sviluppo aggiunto dedotta da una oculata scelta di 
funzioni test sia facilmente verificabile la convergenza del procedimento di omogeneizzazione. 
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9 Omogeneizzazione 


Adesso siamo pronti ad utilizzare tutti i vari risultati accumulati per provare la convergenza 
del procedimento omogeneizzazione. 

Dunque, assegnata a(z) € Sxxp, consideriamo il seguente problema di Dirichlet per i valori 
al contorno sul sottoinsieme aperto e limitato G di RN+! 


[g=f ind, 
{ u € Ho(6), (165) 


dove f € L?(G). 
Consideriamo poi il seguente problema omogeneizzato 


Lhmu=f inG, 
{ ue HE(G), FRE 
dove con £L*°" intendiamo |’ operatore di Kolmogorov-Fokker-Planck 
9 
pere 5, ds.(af?"d,,) +(r,BD)-d (167) 
i,.j=1 
i cui coefficienti glo, per î,j = 1,...,9, sono quelli relativi alla matrice A*°" indicata nella 


Proposizione 15 oppure, il che è lo stesso, quelli forniti dalla relazione (118). 
Naturalmente è necessario chiarire in che senso intendiamo che una funzione u € #5°(G) sia 
soluzione debole del problema (166). Allora diciamo che w € H5°(G) è soluzione debole di 
(166) se verifica la seguente identità integrale | 


fate Desa), Ds, (A)dz+ | u(2)(Y%(2))dz = -/ ÎS(2)0(2)dz (168) 
G G S 


per ogni w(z) € C$°(G). Facciamo presente che al momento non sono disponibili risultati 
di esistenza e di unicità per la soluzione debole del problema (166) come appena definita. 
Possiamo adesso enunciare il seguente Teorema che contiene il risultato finale e conclusivo 
del nostro lavoro. 


Teorema 16 Sia u € Ho(G) la soluzione debole, per e fissato, del problema (165). Allora, 
a meno di una sottosuccessione, 


u + u in L?(G)-debole (169) 
dove u € HE°(G) è una soluzione debole del problema (166). 


Prima di passare alla dimostrazione notiamo che l’ esistenza di una soluzione debole del 
problema omogeneizzato è garantita, nel nostro caso, dalla esistenza di una successione ad 
essa convergente. 

Dimostrazione. Cominciamo.col far presente che grazie alle relazioni di convergenza (136) e 
(152) ed alle proprietà del pavage intrinseco è immediato verificare che 


87(2, D(e')z) + di(2, D(e7!)z) in L?(G)- debole, (170) 
dove, per j = 1,2, 97 sono le funzioni definite rispettivamente dalle (139) e (155) mentre 


ile, DEN) = DANTE 


Ty 
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A sl TE AO 
2(2, 1)2) = 4 1)z o 3 
BEDA LRD 


Allora, utilizzando la notazione distribuzionale e scegliendo come famiglia di funzioni test la 
vs; definita nella (133), il problema (127) implica il seguente 


(uf, (Lou) vi) = (f, un) vie Ho(G). (171) 


Adesso avendo presente lo sviluppo aggiunto (164) possiamo riscrivere la relazione precedente 
nella forma 


(fi v?) = (us, (Liu) "vs) = 
= (Un: (£5°7) 00) +e(u5, 97) +e*(ut, 92) +e°(u5:.93) +e%(ut, 91). 


Grazie alla regolarità di cui godono le funzioni 97 (2,6) sia nella variabile macroscopica che 
in quella microscopica e ricordato che 


us Ilex < e 
dove c costante è una costante indipendente dai parametri 7 ed e, si ottiene facilmente che 
Kn: 97) SG; 


per j = 1,...,4 dove le costanti c; sono indipendenti tanto dal parametro 7 quanto dal 
parametro e. Pertanto, grazie a queste considerazioni, si ha che 


Iuns (Liub) 05) — (us, (CAP) vo) < Ce, (172) 
dove C' è una costante che non dipende né da n né da e. Quindi 
IS, vr) — (uf, (CAM)*vo)] < Ce; (173) 


grazie alla (170), alla (157) ed alla Proposizione 14, passando al limite nella precedente per 
np+ 0, si ha 
KA, v°) Di (u*, (£L*°9)"vo)] £ Ce, i (174) 


e poiché 
0 + vo in L?(G)-debole 


mentre, a meno di una sottosuccessione, 
uu in HE (G)- debole, 
passando al limite nella (174) per € + O si ottiene 
(u, (L*")° vo) = (f, vo) (175) 


ovvero u € H5°(G) è una soluzione debole del problema omogeneizzato. 
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